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ABSTRACT 

The young high-eccentricity binary DQ Tau exhibits powerful recurring 
millimeter-band (mm) flaring attributed to collisions between the two stellar 
magnetospheres near periastron, when the stars are separated by only ~ 8 R*. 
These magnetospheric interactions are expected to have scales and magnetic 
field strengths comparable to those of large X-ray flares from single pre-main- 
sequence (PMS) stars observed in the Chandra Orion Ultradeep Project (COUP). 
To search for X-rays arising from processes associated with colliding magneto- 
spheres, we performed simultaneous X-ray and mm observations of DQ Tau near 
periastron phase. We report here several results. 1) As anticipated, DQ Tau was 
caught in a flare state in both mm and X-rays. A single long X-ray flare spanned 
the entire 16.5 hour Chandra exposure. 2) The inferred morphology, duration, 
and plasma temperature of the X-ray flare are typical of those of large flares from 
COUP stars. 3) However, our study provides three lines of evidence that this 
X-ray flare likely arises from colliding magnetospheres: the chance of capturing 
a large COUP-like flare within the span of our observation is small; the relative 
timing of the X-ray and mm flares indicates the Neupert effect and is consis- 
tent with a common coronal structure; the size of the emitting coronal structure 
(4 — 5 R*) inferred from our analysis (which is admittedly model-dependent and 
should be considered with caution) is comparable to half the binary separation. 
4) The peak flare X-ray luminosity is in agreement with an estimate of the power 
dissipated by magnetic reconnection within the framework of a simple model of 
interacting magnetospheres. 
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Introduction 



T-Tauri stars generally show highly el evated levels of X-ray activity arising mostly from 



violent magnetic reconnection events (e.g. iFeigelson fc Montmerldll999l ). This strong X-ray 



emission has far-reaching implications for the physical processes in the circumstellar environ- 
men t, the formation of planetary systems, and t he evolution of protoplanetary atmospheres 



e.g. iGlassgold et al.ll2005l : IFeigelson et al.ll2009h . Recent X-ray surveys of nearby pre- main 



sequence (PMS) stellar populations give detailed insight s into T-Tauri magne tic flaring; these 



include the Chandra Orion Ultradeep Project (COUP; IGetman et al.ll2005l ) and the XMM 



Newton Extended Survey of Taurus (XEST; iGiidel et al.l l2007aj ) . Astrophysical studies of 



the properties of individual flares and statistical studies of many flares, from both the COUP 
and XEST observations, reveal that most events are similar to solar magnetic flaring, but 



peratures up to 50 times higher (e.g. 


7 avata et al. 


2005: 


Woll 


t et al. 


2005; 


Flaccomio et al. 


2005; 


Stassun et al.ll2006 


Maggio et al. 


2007; 


Caramazza et al. 


2007 


;] 


Jranciosini et al.ll2007 


)• 



Accretion shocks may contribute only a small fraction t o the total X-ray emission from T- 



Tauri stars, in the form of soft X-ray excess emission (e.g. iTelleschi et al.ll2007bl ; IGiidel et al. 



2007bh . 



COUP, the deepest and longest (13.2 continuous days) X-ray observation of a young 
stellar cluster, provided a unique opportunity to stu dy relatively rare (t ypically 1 flare per 
week per star) big X-ray flares from T-Tauri stars (jFavata et al.l 120051 1 . More recently, a 
detailed and systematic st udy of > 200 b i g flares from > 150 bright PMS stars detected in 
the COUP observation by IGetman et al.l ( l2008ai lbl) [hereafter G08a and G08b] shows that 
they are the most powerful, longest, and hottest stellar flares corresponding to the largest 
known stellar X-ray coronal structures, reaching up to several stellar radii in both disk- 
bearing (Class II) and diskless (Class III) systems. The associated large-scale magnetic 
fields (assuming a dipolar geometry) have an equipartition strength B = 0.05 — 0.3 kG in 
the outer loop region, consistent with optical Zeeman measu rements of surface field strengths 
around 2 — 3 kG in magnetically active T-Tauri stars (e.g. I Johns-Krulll l2007t iDonati et al. 



20081 ). G08ab also provide observational evidence for magnetospheric truncation by a disk in 
a Class II system, and for the ability of X-ray loops to withstand centrifugal forces in rapidly 
rotating Class III systems. G08ab propose the COUP sample of flares as possible enhanced 
analogues of very rare Solar Long Decay events (LDEs) associated with X-ray arches and 



-3 - 



streamers. 



While observations of PMS stars at millimeter wavelengths are typically used to study 
steady thermal emission from dust in their protoplanetary disks, transient gyrosynchrotron 
and syn chrotron continua from flares with spectral p eaks in the GHz-THz range can also 
be seen (IPriest fc Forbes! 120021 ; Kaufmann et al.lll986). Long-term radio variability in older 
Class III PMS stars has been known for some t ime (IGarav et al.lll987l) and generally does 
not show correlations with X-ray variability (e.g. lForbrich et al.ll2007l ). Short-term (hours- or 
day-long) radio outbursts are occasionally seen from PMS systems, for exampl e: a remark- 



ably powerful millimeter flare from a heavily absorbed Class III s ystem GMR-A (jBower et al. 



20031 ) ; recurring flares from the Class HI binary V773 Tau A (IMassi et al.l 120081 ) and the 



Class II binary system DQ Tau (ISalter et al 



stars in the Corona Australis cloud (I Choi et al 



2008 



2010( 1: IRS 5b and IRS 7A Class I proto- 



2009h: and a poorly c haracterized system in 



Orion's embedded BN/KL star forming region (IForbrich et al.l 120081 ) . Prior to the current 
study, only one of these cases, GMR-A, was simultaneously observed with a modern X-ray 
telescope, by a coincidence with the Chandra-COXJP observation. The GMR-A mm flare was 
associated, though not exactly simulta n eous, with several days of complex big X-ray flares 
feower et all I2OO3I : iFuruva et aD I2OO3I : iFavata et aJ l2005h . The GMR-A star is believed 
to be a single star with a strong coronal magnetic field, while V773 Tau A and DQ Tau 
are close binary (or triple) systems with component separations at periastron of 30 R* and 
8 R*, respectively. For these multiple systems the magnetic reconnection events have been 
attributed to interacting magnetospheres. 

The DQ Tau binary is especially useful in regards to possible X-ray emission. This is 
a non-eclipsing, double-lined spectroscopic binary, comprised of two relatively equal-mass 
(equal-radius) PMS stars of M ~ 0.65 M (R ~ 1.6 R ) with spectral types in the range of 
K7 to Ml, a rotational period of P ~ 3 days for bo th stars, and a robust orbital period of 
15.804 days (IMathieu et al.lll997l ; iBasri et al.lll997l ). Its highly eccentric orbit (e = 0.556) 
exhibits a periastron separation of only ~ 8 R* (~ 13 R©). The spectral energy distribution 
(SED) of DQ Tau is fa irly typical of a Clas s II system, fit by a large circumbinary disk of 
about 0.002 — 0.02 M (IMathieu et al.l 119971 ) . For more than 65% of periastron encounters, 
the system experiences optical brighten ings as a result of variable and irregular accretion 
(IMathieu et al.l Il997t IBasri et al.l Il997l ). The periastron separation is expected to induce 
magnetospheric interactions at scales and magnetic field strengths comparable to those in- 
ferred for the COUP sample of big flares (G08b). We thus performed simultaneous X-ray 
and mm observations of the orbital segment of DQ Tau around the peaks of the previously 



detected mm flares (ISalter et al.ll2008l ) to search for X-ray emission arising from processes 
associated with colliding magnetospheres. 
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We report here a Chandra detection of a long X-ray flare accompanied by mm activity 
cl ose to the per i astron passage of DQ Tau. The mm observations are discussed in more detail 



m iSalter et al.l (120101 ) . The Chandra data are described in §2.1( and the treatment of mild 
photon pile-up in the observation is given in §2.21 Archived XMM-Newton X-ray observations 
of DQ Tau at an orbital phase away from periastron are presented in §2.31 Time-integrated 
Chandra and XMM — Newton spectra are compared in §3.21 The Chandra flare analyses 
and the derived flare loop length and loop thickness (within the framework of a single- 
loop model) are presented in §3.31 and §3.41 A comparison of the Chandra flare with the 
coincident mm flare is provided in §HJ A comparison with the COUP sample of big flares is 
given in £j5j We end in §|6] with a discussion of the applicability of the single-loop approach 
(including information from Appendix |A]on an X-ray analysis of the multi-loop solar X-class 
flares), our observational findings, and their implications for the origin of the X-ray emission, 
energetics, and loop geometry. 



2. X-ray Observations and Data Extraction 



2.1. Chandra Data 



Th e observation of DO Tau was obtained on January 11 — 12, 2010 with the ACIS 
camera ( IGarmire et al.ll2003l ) on-board Chandra ( IWeisskopf et al.ll2002[ ) as a combined Guest 
Observer and Guaranteed Time Observation (Obsld No. 10992, Co-P.I.s: K. Getman & G. 
Garmire). The observation was set up based on the expectation of an X-ray flare close to 
the orbital phase of $ = 0.98, whi ch corresponds to the peak of the April 2008 large sub- 
mm flare ( palter et al.l 120081 . 120101 ) . The total exposure time of the Chandra observation 
is ~ 59 ks, with no data losses or background flaring due to solar activity. The start and 
end times of the observation are January 11, 2010 at 14:02:43 UT (JD 2455208.09) and 
January 12, 2010 at 06:31:17 UT (JD 2455208.77), respectively. This time period covers 
the orbital segment of $ = 0.95 — 0.99 assuming the orbital parameters, such as time of 
periastron pas s age JDp = 2449582.54 and orbital period P = 15.8043 days, determined by 
Mathieu et al.l (119971 ) . To mitigate photon pile-up effects during an anticipated X-ray flare, 
the observation was performed with a 1/8 sub-array of a single ACIS-I3 chip. The aim point 
of the observation was 04 h 46 m 54!2, +16°59 / 35 / /3 (J2000), and the satellite roll angle was 
284? 1. 



Data r eduction follows proce dures similar to those described in detail by iBroos et al. 



(I2010l ) and iTownsley et al.l ( 120031 . Appendix B). Briefly, using the tool acis_process-events 
from the CIAO 4.2 software package, the latest calibration information (CALDB 4.2.0) on 
time-dependent gain and a custom bad pixel mask are applied, background event candidates 
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are identified, and the data are corrected for CCD charge transfer inefficiency (CTI). Using 
the acis_detect_afterglow tool, additional afterglow events not detected with the standard 
Chandra X-ray Center (CXC) pipeline are flagged. The event list is cleaned by "grade" 
(only ASCA grades 0,2,3,4,6 are accepted), "status", "good-time interval", and energy fil- 
ters. The slight point-spread function (PSF) broadening from the CXC software position 
randomizations is removed. 



Using the ACIS Extract (AE) software package (jBroos et al. 2010 ). DQ Tau photons 



are extracted within a polygonal contour enclosing ~ 99% of the local point spread function 
(PSF). The background is measured locally in a source- free region (Figure Hk). More than 
6200 DQ Tau net counts in the full (0.5-8 keV) band are detected. The AE package was also 
used to construct source and background spectra, compute redistribution matrix files (RMFs) 
and auxiliary response files (ARFs), construct light curves and time-energy diagrams, per- 
form a Kolmogorov-Smirnov (K-S) variability test, compute photometric properties, correct 
spectra and light curves for light pile-up ( §2.2j) . and perform automated spectral grouping 
and fitting of time-resolved data ( §3.3p . 

With a fast rise phase (e-folding timescale r r i se = 26.3 ks) and a slow decay phase 
{jdecay = 40.9 ks), a single long X-ray flare spanned the entire ~ 16.5 hour Chandra exposure 
(Figure [lb). During the observing period the count rate in the full-band (and hard-band) 
increased by a factor of 2 (and 3), peaking at about t = 27 ks ($ ~ 0.97) after the start of the 
Chandra observation. The evolution of the median energy of X-ray events in the full-band 
(adaptively smoothed) shows a slower rise and a faster decay than the light curve, reaching 
a ~ 5 ks wide maximum-value plateau at t ~ 23 ks, i.e. 4 — 5 ks before the count rate peak 
(Figure [Tb,c). The X-ray count rate and median energy serve a s photometric surrogates for 



emission measure and plasma temperature (IGetman et al.ll2010l ). There is thus an indication 



for a time delay between the temperature peak and the emission measure peak (see §3.31 for 
an explanation of the effect). 



2.2. Chandra Data Pile-up Analysis 

Photon pile-up occurs when two or more photons are incident on a single detector region 
(for Chandr a- ACIS this is typically a 3 x 3 pixel region) during a single CCD frame. Two 
or more photons are thus detected as a single event. Some of these pile-up events mimic 
cosmic rays and are rejected by the on-board processing, while others are telemetered as 



1 The A CIS Extract software package and User's Guide are available at 
http : //www. astro .psu. edu/xray/acis/acis_analysis .html. 
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valid events but with spuriously high energies^. Despite the use of the 1/8 sub-array ACIS 
mode with the reduced CCD frame of 0.5 s (6.4 times shorter than the nominal frame), the 
Chandra data of DQ Tau suffer minor pile-up. Assuming a single-temperature plasma with 
kT ~ 2 keV and an average X-ray column density of N H ~ 2 x 10 21 cm~ 2 ( §33} , PIMMjfl 
gives an estimate on the pile-up fraction (the ratio of the number of frames with two or 
more events to the number of frames with one or more events) of only 3% at the peak of 
the Chandra flare. Nevertheless, below we adopt two independent approaches to correct 
the data for pile - up: m odeling pileup within the context of spectral fitting using the pile-up 



model of iDavisI ( 120011 ). and reconstructing a pile- up free spectrum using an experimental 
Monte Carlo approach of Broos et al. (2010, in prep.). Both methods are applied to the 
forty five 1000-count overlapping time segments defined in §3.31 



In the first approach, for each of the 45 time segments, the pile-up model of Davis is 
applied along with fixed two-temperature soft component and thawed one-temperature hot 
(flare) component plasma models subject to constant absorption, as defined and discussed in 
§3. H and §3.31 According to "The Chandra ABC Guide to Pileup", the following parameters 
are left frozen at their default values: the maximum number of photons to pile-up in a single 
frame max p h = 5, the grade correction for single photon detection g = 1.0, and the number 
of independent 3x3 pixel pile- up islands in the source extraction region nregions = 1. 
The frame time parameter frtime is set to that of the Chandra observation's frame time 
of EXPTIME = 0.5 s divided by the fractional exposure FRACEXPO = 0.995. The 
grade migration parameter (a), and the fraction of events in the source extraction region 
to which pile-up is applied (psffrac), are varied within the ranges of a = [0.5 — 1.0] and 
psffrac = [0.90 — 0.95], respectively. 

In the second approach, an experimental Monte Carlo forward-modeling method to re- 
construct unpiled DQ Tau spectra from piled-up ACIS data is applied. The essential features 
of this pile-up reconstruction method are briefly described here. A thorough description of 
the method is given by Broos et al. (2010, in prep). A non-physical "nuisance model" with 
many free parameters feeds an input spectrum to the MARX mirror-detector simulator^ 
which produces photons that have the correct Chandr a- ACIS PSF. A physical model of the 
ACIS CCD produces a piled simulated spectrum. The nuisance model is iteratively adjusted 
until the piled simulated spectrum is similar to the observed spectrum. The simulation is 



2 Detailed information on a Chandra- ACIS pile-up is given, for example in "The Chandra ABC Guide 
to Pileup" , http : / / cxc .harvard. edu/ciao/download/doc/pileup_abc . pdf 

3 http : //asc . harvard . edu/ toolkit /pimms . j sp 



http : //space . mit . ed u/CXC/MARX/ 
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then run one more time, with pile- up disabled (i.e. exactly one photon arrives per frame). 
The resulting event list is claimed to be similar to what ACIS itself would have produced if 
pile- up were not present, and thus one could attempt to fit it with a physical model. For 
each of the 45 time segments, reconstructed unpiled spectra of DQ Tau are then fit with fixed 
two-temperature soft and thawed one-temperature hot plasma component models subject to 
constant absorption, as defined and discussed in §3.11 and §3.31 

For all 45 time segments, flare plasma temperature and emission measure inferred by 
the two methods are consistent within la. Both methods indicate very light pile-up with 
< 4% difference between the observed count rate and that expected if pile-up effects were 
not present (filled versus open circles in Figured]). Further in the flare analysis (§3J) we adopt 
the results from the Broos et al. method. 



2.3. XMM-Newton Data 



We use archived, publicly available XMM-Newton (IJansen et al.ll200ll ) observations of 
DQ Tau. As part of the program aimed at studying X-ray properties of EX Lupi-type out- 
burst (EXOR) stars, the EXOR star DR Tau was observed by XMM-Newton on February 13, 
2007 (Obsld No. 0406570701, P.I. G. Stringfellow) . DQ Tau is 3 3' north-west of PR Tau 
within the field of view of the XMM-Newton EPIC Jstriider et al.lboOll ; iTurner et aJboOlh 
detectors (Figure [2^). Both EPIC-MOS and EPIC-PN cameras were configured in full win- 
dow mode with medium and thin filters for MOS 1,2 and PN cameras, respectively. The 
total exposure time for both EPIC-MOS 1,2 detectors is 12.6 ks, and it is 11 ks for the EPIC- 
PN. The observation start and end times, 04:20:12 UT (JD 2454144.68) and 07:50:28 UT 
(JD 2454144.83) respectively, correspond to a time period covering the $ = 0.66 — 0.67 
orbital segment of DQ Tau. 

Our data reduction follows procedures described in the Common and EPIC related 
Science Analysis Software (SAS) Threads^. Briefly, the data were reprocessed with the 
SAS version 9.0.0. The emproc and epproc meta-tasks were run to obtain calibrated and 
concatenated event lists for EPIC-MOS and EPIC-PN detectors, respectively. Short intervals 
of flaring particle background of 100 s (400 s) long around the observation time of t = 6 ks 
for MOS (PN) detectors were identified and removed. The event lists were further cleaned 
with the PATTERN, FLAG, and energy filter fl 



5 The SAS Threads are available at http : //xmm. esac . esa. int/sas/ current/documentation/threads/ 
6 The following niters recommended in http : //xmm . vilspa . esa . es/docs/documents/CAL-TN-0018 . pdf 
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The DQ Tau photons were extracted within a 33.5" radius circular region (~ 90% 
encircled PSF). For the MOS 1,2 images, background was measured on the same CCD away 
from the source; for the PN image, a background region was chosen from a source-free region 
on a neighboring CCD at the same distance to a corresponding readout node (Figure^). 
More than 4100 DQ Tau EPIC net counts in the (0.2 — 10) keV band were detected. The 
SAS tasks evselect and epatplot suggest that the EPIC image of DQ Tau is not affected by 
pile-up. 

Following the SAS threads, the SAS tasks evselect, backscale, rmfgen, arfgen, epiclccorr 
and a number of tools from HEASOFT version 6.8 were used to construct source and back- 
ground spectra, compute RMFs and ARFs, and construct light curves. During the XMM 
observation, DQ Tau showed no significant variations in its X-ray light curve (Figure [2b). 



3. Chandra Flare Analysis 
3.1. Time-integrated Chandra spectrum 



The time- integrated Cha/nd/ra spectru m is modeled here with a three-temperature APEC 
plasma emission model ( jSmith et al.ll200lf ). The APEC model is also employed in the time- 
resolved flare spectroscopy ( §3.3)) . However, for a better modeling of individual spectral 
features, comparison bet ween the Chandra and XMM-Newton data in §3.21 is performed 
using the VAPEC model JSmith et al-lboOlh . 



The fit to the time- integrated Chandra spectrum is performed using x 2 statistics. The 
light Chandra data pile-up ( §2.2)) is ignored here. In the fit, elemental abundances are frozen 
at 0.3 ti mes solar, which has been previous l y suggested as typica l values for young stellar 
objects (llmanishi. Kovama. fc Tsuboil 1200 ll ; iFeigelson et al.l 120021 ). Solar abundances are 
taken from Anders fc GrevesseT(1989|). X- ray absorption is modeled using the atomic cross 



fitting package version 12.5. In (jArnaudlll996l ). 



sections of Morrison &: McCammon (jl983 ) . T he fit is performed with the XSPEC spectral 



The X-ray spectral fitting of the DQ-Tau data is ambiguous. At least two qualitatively 
different spectral model families give similarly good values of reduced xl ~ 0.S@ (see Fig- 



are used. For both the EPIC-MOS and EPIC-PN detectors, the energy filter of (0.2 - 10) keV is applied. For 
the EPIC-PN detector, excluding border pixels (FLAG == 0), single- and double-pixel events (PATTERN 
<=4) flagged as "good" (#XMMEA_EP) are retained. For the EPIC-MOS detector single-, double-, triple-, 
and quadruple-pixel events (PATTERN <=12) flagged as "good" (#XMMEA_EM) are retained. 

7 Ambiguity in spectral fitting of hundreds of bright X-ray young stellar objects was previously reported 
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ure[3^,b with analogous VAPEC fits). The best fit parameters for the first model family are 
around kT x = 0.7 keV, kT 2 = 1.9 keV, kT 3 = 4.3 keV, N H = 1.3 x 10 21 cm" 2 (c orresponding 
to a v isual extinction of Ay ~ 0.8 mag using the gas-to-dust relationship of IVuong et al. 
(120031 )). and for the second family are around kTi = 0.3 keV, kT 2 = 0.9 keV, kT 3 = 3.5 keV, 
Nh = 3.2 x 10 21 cm -2 (Ay ~ 2 mag). The second solution gives softer temperatures and 
higher absorption. The high temperature model component (kT 3 ) is assumed here to be 
fully associated with the X-ray emission from the flare. A summary of the low temperature 
model components (kTi and kT 2 ) for both solutions is given in Table [TJ Both seem phys- 
ically reasonable. For example, in the first solution, the temperatures (kT\ and kT 2 ) and 
the ratio of emission measures (E M 2 / EMy ) are consistent with those for COUP PMS stars 
(Table 1 in iGetman et al.l 120101 ) . In the sec ond solution, fcXi ~ .3 keV might represent 



a soft excess emission due to accretion (e.g. ITelleschi et ajj|2007bj) ; DQ Tau is known to 



experience pulsed accretion flows near periastron (IBasri et al.lll997l ). For both solutions, the 
absorption-corrected total-band (0.5-8 keV) luminosities of the low temperature component 
(kTi and kT 2 ) lie within the loc us of XEST Class II stars on the X-ray luminosity versus mass 
diagram (ITelleschi et al.ll2007ai ). Finally, the X-ray column densities from both solutions are 
consistent with t he range of DQ Tau v isual extinction esti mates given in the literature, from 
A v = 0.5 mag JCohen fc Kuhilll979l ) to Ay = 2.1 mag jstrom et alill989l ). Thus, in the 
following analysis of the Chandra flare we will consider both solutions, which we will refer 
to as low temperature plasma component Model 1 and Model 2. 



3.2. X-ray Emission at a Different Epoch and Phase 



Compared to the emission near periastron ($ = 0.95 — 0.99) observed in January 2010 
by Chandra, the emission away from periastron ($ = 0.66 — 0.67) observed in February 2007 
by XMM-Newton is rather constant (Figure [2b) and, as shown below, softer. The time- 
integrated Chandr a- ACIS and XMM- Newton-FjPld^ spectra are compared in Figure O Both 
spectra are fit with a two- or three-temperature VAPEC plasma emission model with indi- 
vidual elemental abundances set to values typical for PMS and for extremely active zero-age 
main-sequence (MS) stars, as specified in Giidel et al. ( 2007 a)^l. The XMM-Newton spec- 



bv lGetman et al.1 (|2005l ) 



8 This is a combined PN+MOSl,2 spectrum created following the SAS thread 
http : //xmm. esac . esa. int/sas/ current /document at ion/threads/ epic_merging . shtml 

9 With respect to the solar photospheric abundances of Anders fc Grevessel ( 19891 ). individual elemental 
abundances adopted in our spectral fits are: C=0.45, N=0.788, 0=0.426, Ne=0.832, Mg=0.263, Al=0.5, 
Si=0.309, S=0.417, Ar=0.55, Ca=0.195, Fe=0.195, Ni=0.195. 
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trum can be successfully fit by the low temperature plasma component Model 1 with no need 
for a hotter component, but with a flux 1.5 times higher than that of the low temperature 
plasma emission in the Chandra spectrum (panels (c) versus (a) in Figure [3]) • Alternatively, 
the XMM-Newton spectrum can be fit with the low temperature plasma component Model 2 
with a flux 1.7 times lower than that of the low temperature emission in the Chandra spec- 
trum, and an additional moderately hot component (kT 3 = 1.9 keV) with a flux 4.5 times 
lower than that of the Chandra hot component (kT 3 = 3.5 keV) (panels (d) versus (b)). 
Differences in plasma emission between the two epochs, spaced roughly 3 years apart, could 
be simply attributed to evolution of magnetic active regions on the stellar surface. The 1.5- 
fold flux difference between the XMM-Newton spectrum and the Model 1 low temperature 
component emission in the Chandra spectrum has an alternate, less likely, explanation - 
the known suppression of time-integrated X-ray emiss ion in accreting versus non-accreting 
stars (e.g. iFlaccomio et al.ll2003l ; iPreibisch et al.ll2005l ); phase- variable mass flows from the 
circumbinary disk with accretion rates 10 — 200 times higher near periast ron than those 
away from periastron are proposed to take place in the DQ Tau system ( jMathieu et al. 



19971 : iBasri et al.lll997fl . 



3.3. Time- Resolved Spectroscopy 

As in G08a, flare spectral modeling is performed here on multiple adjacent overlapping 
time segments. Each time segment is specified by moving a rectangle ("boxcar") kernel of 
variable width through the Chandra event time series to encompass 1000 X-ray counts. Start 
time of each segment is specified to be offset from the start time of the previous segment 
by 0.1 times the duration of the previous segment. Forty five segments are defined. The 
typical width of the smoothing kernel is 7 ks at the peak of the flare and 10 ks at the base of 
the flare (Table [2]). Among the 45 overlapping segments, sets of 5-6 independent segments 
can be identified. Examining dozens of overlapping time segments along the decay phase 
of a flare light curve, rather than just a handful of independent intervals, allow higher time 
resolution which often results in the discovery of a more detailed, often more complex flare 
behavior (Appendix A in G08a). 

Both piled and unpiled spectra ( §2.2|) from each segment are fit in XSPEC with thawed 
one-temperature APEC hot component (flare) and fixed two-temperature APEC cooler com- 
ponent (Tabled]) models subject to constant WABS absorption. As in §3.11 elemental abun- 
dances are set to 0.3 times solar, and fits are performed using \ 2 statistics. Examples of the 
spectral modeling are shown in Figure HI and the unpiled spectral results (plasma temper- 
ature, emission measure, corrected-for-absorption X-ray luminosities) for both Model 1 and 
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Model 2 scenarios are reported in Table [2j 



For Model 1 and Model 2, respectively, Figures and E] give the inferred temporal evo- 
lution of the plasma temperature and the emission measure, as well as flare evolution in 
the logT — log V EM plane. Without pile-up correction (gray points), flare temperatures 
would have been overestimated by 20 — 25% and emission measures would have been un- 
derestimated by 10%. During the first 15-20 ks, the temperature profile shows a plateau 
with a relatively constant plasma temperature around 50 — 60 MK, followed by a drop to 
40 — 50 MK. This might be an indication of a precursor flare event. Afterwards, the tem- 
perature profile is typical of many big flares from the COUP sample of flares: starting with 
a 5 ks period where the temperature rapidly rises to 60 — 90 MK, then in the next 5 ks it 
rapidly drops to 30 — 40 MK, and decreases more gradually to 25 — 30 MK over the following 
> 10 ks. 

A 5 — 8 ks time delay is seen between the temp erature peak and the emission measure 
peak. Such delays are often observed in sol ar fe.g.. iMilkev et al.l Il97ll ; iGuedel et al.l Il996 



and Appendix IA1 of this work), stellar (e.g., iReald 120071 . and references therein), and flares 
from PMS stars (e.g., Table 1 in G08a). The delay can be interpreted as an efficient cooling 
of the plasma by conduction and radiation while th e emission measure is still rising due to 
ongoing evaporation in the loop (IGuedel et al.lll996l ). The effect is applicable to both flares 
from single and multiple loops (see Appendix |A|) . 

The derived slope of the trajectory in the DQ Tau flare's temperature-density diagram 
£ = 1.1 — 1.5 indicates mild or no sustained heating in the framework of a single- loop flaring 
model (see §3.4.11 for more details). 



3.4. Single-Loop Modeling 

3.4-1- Flare Loop Length 



In order to derive the size of the coronal structure associated with the Chandra flare in 



DQ Tau, we employ the time-dependent hydrodynamic model of iReale et al.l (119971 ) [hereafter 
R97] for a single dominant coronal magnetic loop. Those authors establish a formula for 
estimating a loop's half-length L (throughout the text the terms "size" or "length" will be 
used to indicate the loop's half-length) accounting for the possibility of prolonged heating 
during the decay phase. They find 



'Tdecay yjTpk 

3.7 x 10" 4 F(C) 



(1) 
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where L is the half-length of the loop (cm), Tdecay is the flare decay e- folding timescale (sec), 
and T pk is the loop apex flare temperature (K) at the time of the maximum emission measure. 
F(Q is a correction factor for prolonged heating that is a function of the slope ( of the 
trajectory in the temperature-density diagram. In practice, F(() and T pk must be calibrated 
for each X-ray observatory; the slope ( is usually measured in the logT — \og(EM 1 / 2 ) plane 
where EM is the evolving emission measure and EM 1 / 2 is used as a proxy for the plasma 
de nsity. We reproduce the calibration formulas for F(() and T pk derived for Chandr a- ACIS 
bv lFavata et~ZI f|2005h 



and 



F(0 



0.63 



C-0.32 



+ 1.41 



(2) 



0.068 x TV 2 



l pk = U.UUO A l obg (3) 

where T pk is a temperature at the loop apex (K), and T obs is an observed "average" loop 
temperature (K) obtained from our Chandra- ACIS data. The ^ ~ 1.5 corresponds to a freely 
decaying loop with no sustained heating, while ( ~ 0.32 corresponds to a loop with a long 
sustained heating. 

This model is simplistic in a number of ways. It assumes that the plasma has a uni- 
form density with a unity filling factor confined within a single semicircular loop of uniform 
cross-section. Furthermore, the model assumes that this geometry remains unaltered during 
the flare, that energy is efficiently transported along the magnetic field lines of the loop, and 
that there is continuous energy balance between the loop heating and the thermal conduc- 
tion and radiative losses. Despite these limitations, the R97 model has been applied to a 
variety of solar and stellar flares, including flar es from PMS stars i n the ONC and Tauru s 
molecular cloud s t ar-for ming; regions studied by iFavata et al.l (120051 ) ; iGetman et al.l (j2008al ) ; 



Franciosini et al. 



(2007 1). The model has the adva ntage over earlier and simpler cooling loop 



models (e.g. iRosner et al.lll978t ISerio et al.lll99ll ). which neglect reheating during the decay 
phase and thereby tend to overestimate loop sizes. Note, however, the limitations of the 



single- loop approach and its application to the DQ Tau flare discussed in §6.11 

The decay timescale inferred for the DQ Tau flare is Tdecay = 40.9 ks (§2]). In the 
case of the Model 1 (Model 2) the observed "average" loop temperature at the time of 
the peak of the emission measure (time segment # 27 in Table [2]) is T obs = 68.6tg 4 9 ° MK 
(T b s = 48. MK), and the slope of the trajectory in the temperature-density diagram is 
C = 1.5 ± 0.6 (( = 1.1 ± 0.5; see Figures [ME]) . Applying formulas (II]) - (EJ) we find the size 
of the DQ Tau flaring coronal structure L = (5.8 ± 0.7) x 10 11 cm (L = 5.2 ± 0.6 R*) and 
L = (4.4 ± 0.6) x 10 11 cm (L = 4.0 ± 0.6 R*) for the Model 1 and Model 2, respectively. 
This size is comparable to independen tly inferred height s for the mm emitting regions of 
3.7 — 6.8 R* above the stellar surface (ISalter et al.ll2010i ). The loop size of L = 4 — 5 R* 
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i s also comparable t o half the > 8 component separation near the periastron passage 
jMathieu et al.lll997h . 

The R97 model assumes that loop growth is unrestrained by gravity. This is clearly the 
case for DQ Tau flare, where the derived half-length of the flaring loop is much shorter than 
the pressure scale heighfl H r , both near the stellar surface (H p ~ 7 — 16 R*) and at a few 
stellar radii from the surface (H p > 100 R*). 



3.4-2. Cooling Timescales and Flare Loop Thickness 



The analysis described in §3.4.11 only provides the loop length. The ratio between the 
loop's cross-sectional radius and the loop length (f3 = r/L), which indicates the thickness of 
the loop, remains to be determined. For the cool (T ~ 1 MK) solar active region coronal loops 



spati ally resolved by the Transition Region and Coronal Explorer ( TRA CE, lHandy et al. 



19991 ). j3 ranges from 0.01 to 0.06 for loops with a half-length of L < 10 Mm, and (3 falls 



in th e narrow range 0.007 — 0.02 for larger loops (Table 3 and §4.2 i n lAschwanden et al. 



2000 



ii 



For giant coronal flaring structures found in ONC PMS stars, iFavata et all J2005h 



and G08a assumed /3 = 0.1. This assumption might not hold (section 3.2 in 



Favata et al. 



20051 ). To check the applicability of the Reale method to big PMS flares, IFavata et al. I (boosh 



performed detailed hydrodynamic simulations of the big flare in COUP star # 1343. Their 
loop simulations with f3 = 0.02 were in good agreement with the data. 

Similar to the DQ Tau flare with no (Model 1, ( = 1.5) or mild (Model 2, ( = 1.1) 
sustained heating, the COUP 1343 flare loop is freely decaying with no sustained heating 
(C > 1-5; G08a). Assuming similarity in flare decay cooling processes between the DQ Tau 
and COUP 1343 flares and using the COUP 1343 flare as a testbed, further comparison of 
the cooling timescales with the light curve decay timescales for both flares allows estimation 
of the thickness of the DQ Tau loop. 

The dominant cooling mec hanisms of a coron al plasma are thermal conduction and 
radiation with timescales of (e.g. ICargill et al.lll995l ): 



ln The pressure scale height is H p = (kbT i, s ) / \fj,m p g) . The gravitational acceleration near the stellar surface 
of either DQ Tau component, g, is 0.25 x g & = 6900 cm s 2 . m p is the proton mass and \xm p = 0.6 x m p is 
an average mass of a particle in a fully ionized solar plasma. 

u The loop thickness here is estimated as (3 — w/2/L, where w is the width of a bundle of loop threads 
and not a n individual thread , and L is the loop half-length; w and L measurements are taken directly from 
Table 3 of lAschwanden et alJ (|2000h . 
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3n P khTL 2 



3n P khT 



7~rad 



(4) 



KoT7/2 ,-raa n 2p( T ) 

where fc& is the Boltzmann constant, n e is the electron density, T is the me asured plasm a 



temperature, L is the loop length, k is the coefficient of thermal conductivity (jSpitzer 
optically thin plasma (e.g 



19651 ) 



Rosner et al. 


1978: 


Me we et al. 


1985) 



expressed as functions of independent quantities already derived in our previous analyses, 
such as, the loop length L ( §3.4. lj) . temperature T, emission measure EM, X-ray luminosity 
in the wide 0.01 — 50 keV band £x,o.oi_5o (Table [2]), as well as the quantity of interest here, 
the loop thickness 0. At any given time t: 



2 t 3 



V = 2n(3 z L 



(5) 



n e (t) 



EM{t) 



V 



it) = n e (t) 2 VP(T) 

where V is the loop volume, so that 



(6) 
(7) 



3k h 



EM(t)L 
2^ ' 



T~rad{t) 



3k b T{t)(3 

^X,0.01_5o(^) 



^EM{t)2TxU>. 



For the COUP 1343 flare, the evolution of its T, EM, and £x,o.oi_5o as derived by 
G08a is shown in Figure [7k. Figure [T b-f compares coo l ing timescales obtaine d from equation 
(JSJ) for the loop length of 4.8 R* (IFavata et all 120051 ; iGetman et all l2008ai ) and a range of 
(3 with the observed light curve decay timescale of Td eca y = 24 ks (G08a). For solar and 
stellar flares, if sustained heating is mild or absent, it is anticipated that from the end of 
the heating phase (indicated by the temperature peak) to the time of peak density (the 
emission measure peak) the plasma cooling is governed by conduction; at the density peak 
the radiation c ooling and c onduction times become equal; afterwards the radiation cooling 
dominates (e.g. lRealdl2007l ). During the flare decay the combined conduction and radiation 
cooling time r th (l/r th = l/r con + l/r rad ) is naturally expected to be somewhat close to 
the observed flare decay time r^ecay That is precisely what happens in the evolution of the 
COUP 1343 cooling timescales with a choice of /3 = 0.02 (panel c). In this case, during 
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the flare decay the r t h differs from Tdecay by no more than a factor 1.3, with the conduction 
and radiation times matching at the emission measure peak t = 230 ks and with radiation 
losses dominating afterwards. With a choice of (3 < 0.02, is systematically lower than 
T decay (panel b), while with (3 > 0.02, r th is systematically higher than r decay (panels d-f) 
with radiation losses becoming negligible at (3 ~ 0.1 (panel f). For the COUP 1343 flare 
our cooling timescale analysis thus gives = 0.0 2 and is in complete agreement with the 



results of the detailed simulations of iFavata et al.l (120051 ). This validates our use of cooling 



timescales for the estimation of the loop thickness. 

This method has also been applied to the DQ Tau flare; the results of the modeling 
with the Model 1 are shown in Figure Panels b-e compare cooling timescales obtained 
from equation for a loop length of 5 R* ( §3.4.ip and a range of (3 with the observed light 
curve decay timescale of r decay = 40.9 ks (§2]). With a choice of (3 = 0.03 (Figure Et) the 
evolution of cooling timescales in DQ Tau follows the correct pattern seen in Figure [Tt for the 
COUP 1343 flare. With a choice of (3 < 0.03, r t h is systematically lower than Tdecay (Figure 
[8b), while for (3 > 0.03, r t h is systematically higher than Tdecay (Figure [HH-e). The results of 
the cooling timescale analysis are also in agreement with the loop length analysis (§2]): with 
a choice of loop length L <C 5 R*, r th is systematically lower than Tdecay (Figure [Sf), while 
for L > 5.5 R*, r t h is systematically higher than Tdecay (Figure 02-i). The model L = 4.5 R*, 
(3 = 0.03 is also plausible (Figure [Hg). For the DQ Tau flare our cooling timescale analysis 
thus gives (3 = 0.03 in the case of Model 1 (Figured]), and (3 = 0.06 in the case of Model 2 
(results are not shown). The Model 2 assumption gives a shorter and thicker flaring loop 
with the twice the volume of a Model 1 loop. 



4. Comparison of the Chandra X-ray and IRAM mm Light Curves 



For most of the Chandra observing time the mm flux of DQ Tau was monitored by a 
number o f mm-band i nterfe rometers: IRAM, CARMA, and SMA. The mm data are pre- 
sented by lSalter et al.l (120101 ) . Figure^ compares our X-ray data with the 90 GHz (3.3 mm) 
IRAM data. The IRAM data suggest the presence of at least two microwave flares. The 
end of the captured decay portion of the first microwave flare is close in time to the peak 
of the X-ray flare- Such a light curve behavior is consistent with the Neupert effect (e.g. 
Veronig et al1l2002l ). 



More precisely, the Neupert effect is a correlation between the time- integr ated microwav e 
(or hard X-ray) light curve and the rising portion of the soft X-ray light cur ve (|Neupertlll968l ). 



The effect has been observed in many solar (e.g. iDennis &: Zarrol Il993l ) and some stellar 



flares (e.g. iGiidel et al.l 120021 ). and can be interpreted by the classical non-thermal thick- 



- 16 - 



target model (jBrownlll97ll ; iLin fc Hudsonlll976l ). In this model, electrons are accelerated to 
high energies and spiral downwards along the magnetic field lines of the coronal loop. The 
high-pitch angle (high energy) electrons are trapped in t he coronal re gion of the loop, with 



the highest-energy population emitting microwaves (e.g. lDennislll988l ). The low-pitch angle 
electrons propagate to the chromosphere. When they reach the chromosphere, they produce 
non-thermal hard X-rays and drive chromospheric evaporation, which fills the loop with hot 
plasma emitting in soft (Chandra band) thermal X-rays. 



4.1. Modeling of Thermal Energy Profiles 

Quantitative evaluation of the relation between the IRAM and Chandra data o f DQ Tau 



can be deduced within the framework of the generalized Neupert effect introduced in lGuedel et al. 



(119961 ) . According to this formalism, the temporal profile of the total thermal energy E(t) of 
the plasma in the flaring loop can be expressed as the convolution of the kinetic energy influx 
by non-thermal electrons into the chromosphere (assumed proportional to the observed radio 
flux) with radiative and conductive cooling (assumed to follow an exponential decay) 



E(t) = a f F R (t')e-^ )/fth dt', (9) 

Jt 

where a is an assumed constant conversion factor between the thermal energy flux from 
electrons injected into the chromosphere and the observed radio flux F R , fy, is the time- 
averaged cooling time of the thermal plasma in the loop, to is the start time of the radio 
flare, and t' < t. 

In the case of DQ Tau, the captured late decay phase of the mm flare observed with 
IRAM can be fit to an exponential model 

F R (t) = F R (t 1 )el fL -*V T *, (10) 

where t% = 15.5 ks is the earliest time of the portion of the mm flare decay that was 
observed and r R = 15.9 ks is the decay e-folding timescale during that phase (Figure Ek) • 
The expected thermal energy profile is then 

E(t) = d + aF R {h) f e itl - t ' )/TR e^ t - t ' )/fth dt', (11) 
where the averaged plasma cooling time over the rise phase of the Chandra flare is 
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fth ~ 30 ks (see Figure [Bb) and the C\ is a constant term accounting for the integration over 
the unseen part of the mm flare. The integral in equation (12) can be solved analytically; 
however two constant parameters remain unknown: the expected total thermal energy at 
the end of the unseen part of the mm flare (Ci) and the rate of the injected energy at the 
beginning of the captured decay of the mm flare [a x F^{t\j). 

Equation (12) defines the thermal energy profile using the IRAM data. On the other 
hand, using equations (5) and (6) the thermal energy E(t) = 3n e (t)k b T(t)V = 3k b T(t) v /2nEM(t)(3L 3 / 2 
can be inferred from our X-ray flare analysis; the resulted profile (using X-ray Model 1) is 
shown as the thick solid curve in Figure [9b. Although the energy is subject to significant 
uncertainty, we are not so much interested in the individual values as in the general evolu- 
tion of thermal energy that is observed: an early rise phase (t < 19 ks), when the energy 
is approximately constant around 3 x 10 35 ergs; a late rise phase, when the energy rapidly 
rises to ~ 5 x 10 35 ergs within a 6 ks period; and a flare decay phase, when the energy falls 
to ~ 2 x 10 35 ergs during the remaining 15 ks period. For Model 2 the evolution is similar, 
but the energies are systematically higher by factors of ~ 1.2 — 1.4. 

A family of solutions to equation (12) exists with constant parameter values around 
C\ = 0.5 x 10 35 ergs and a x Fn(ti) = 0.8 x 10 32 ergs s -1 , resulting in thermal energy 
profiles (such as the dashed curve in Figure [9b) that resemble the profile derived from the 
X-ray modeling during the late, but not the early, stage of the X-ray flare rise phased). This 
resemblance of the thermal energy profiles derived independently from the microwave and 
X-ray modeling provides a quantitative indication of a Neupert-like effect. The divergence 
of the dashed and solid curves in Figure [9b for the very early rise phase of the X-ray flare 
(t < 19 ks) could be due to three possible effects that we have not modeled: 1. X-ray emission 
excess from flaring events in other coronal loops; 2. time dependence of the energy conversion 
factor a during the early decay phase of the microwave flare; 3. inaccurate estimation of 
the conduction cooling timescale r con , electron density n e , and thermal energy E during the 
very early stages of chromospheric evaporation. 

Finally, we assume that the second mm flare detected by IRAM is related to a separate 
and independent flaring event in another coronal loop. The weaker X-ray emission associated 
with the secondary event has a negligible impact on the observed decay phase of the primary 
large X-ray flare. 



12 Parameters C\ = 0.5 x 10 35 ergs and a x Fn(ti) = 0.8 x 10 32 ergs s _1 are the values to use in equation 
(12) to match the thermal energy profile from the X-ray modeling with the Model 1. For matching the 
energy profile from the X-ray modeling with Model 2, the best values to use are C\ = 2 x 10 35 ergs and 
a x F R {t x ) = 0.6 x 10 32 ergs s -1 . 
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5. Comparison with COUP PMS and Flares on Older Stars 

Figure [10] compares the DQ Tau X-ray flare properties (indicated with a black box) with 
those of the COUP sample of flares analyze d by G08a (circles), stellar flares from PMS and 



old er active stars com p iled i n iGiidell ( 12004 ) (gray boxes), and solar-stellar trends obtained 
by lAschwanden et all (120081 ) (solid and dotted gray lines). It is important to note that 
the analysis of G08a is limited to COUP sources with > 4000 counts and to flares with a 
peak count rate that is > 4 times the characteristic level. Thus, in Figure [10] the lower 
boundaries of the COUP flare duration, energy, emission measure, and loop size correspond 
to the imposed selection. Panel (a) shows that the rise and decay e-folding timescales of 
the DQ Tau flare lie well within the locus of the COUP sample of flares, indicating similar 
flare morphology. Panels (b) and (e) show that the duration, peak plasma temperature, 
and length of the associated coronal structure of the DQ Tau flare are typical of the COUP 
sample of flares, and are consistent with only the longest and hottest flares from older active 
stars. However, at a given temperature, the DQ Tau flare emission measure is at the lower 
boundary of the COUP and stellar loci (panel (d)). As a result of the relatively low emission 
measure, the DQ Tau flare energy is only comparable to that of the least energetic flares 
from the COUP sample (panel (c)B 

Compared to typical solar flares, the COUP/DQ Tau 
flares are 10 — 100 times longer lasting and hotter, and are associated with 100 — 1000 times 
larger coronal structures. 

The L vs. T b StP k information in panel (e) must be interpreted with care. Since the locus 
of typical solar flares (gray line) represents direct observational measurements of flare length 
scales, and the plotted symbols are model-derived estimates, it is prudent to suspect that the 
relatively large loops inferred for ONC PMS stars (black circles) may represent systematical 
biases introduced by single-loop modeling. However, about half of the 20 flares from older 
active stars (gray boxes) that lie within the locus of solar flares (T b StP k < 40 MK) have been 
modeled using the same single-loop approach of R97; no large upward bias in their inferred 
loop lengths is seen. We caution against concluding that PMS flares follow a different L 
vs. Tobspk correlation than solar flares (grey line, extended), because observational selection 
effects could easily be hiding a very large population of flares that do follow the extended 
solar locus. 

The difference between DQ Tau and typical COUP flare emission measures can be at- 
tributed either to the COUP data selection effect mentioned above, or to the differences in 
flare physical processes. In the latter case, with similar flare morphologies, durations, tem- 



13 In panel (c) the DQ Tau flare energy estimate is given for the (0.5 — 8) keV energy band to be compatible 
with the COUP data. 
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peratures, and loop lengths, the emission measure difference can be explained via distinct 
coronal loop geometries (e.g., a DQ Tau loop that is thinner than typical loops associated 
with flares from the COUP sample) and/or by distinct efficiencies in chromospheric evap- 
oration that in turn can be linked to differences in the magnetic reconnection processes 



6. Discussion 

6.1. Single Loop versus Multiple Loop Scenarios 

Large flares on the Sun often involve arcades of dozens or hundreds of sequentially 
reconnected magnetic loops, which are often obser yed to have different tem peratures with 
the cooler loops lying below the hotter ones (e.g. iReeves &: Warren! 12002k and references 



therein). Despite the morphological complexity of such flares, their X-ray light curves often 
have simple exponential rise and decay shapes0 In view of this solar analogy, the concept 
of the single-loop approach and its application to powerful spatially unresolved stellar flares 
might be questionable. 

R97 have applied their single-loop method t o flare decays of 20 solar M- and C-class 
flares monitored by the Soft X-ray Telescope (SX T, iTsuneta et al.lll99lf ) on board the Yohkoh 
solar observatory satellite ( jOgawara et al.lll99ll ). Loop sizes derived using the method were 
shown to agree with the length of X-ray structures measured from direct inspection of SXT 
images (Figure 6 in R97). In Appendix |A] we apply the approach of R97 to 5 solar X-class 
flares that are clearly associated with arcades of multiple loops seen in EUV TRACE images. 
In contrast to the general view that single loop models tend to overestimate flaring loop sizes 
of complex flare events, for all of our 5 testbed flares the equations of R97 yield a single-loop 
length comparable to or shorter than the lengths of the individual EUV loops measured from 
the TRACE images. 

Here we are unable to give a definite answer to the question of why the loop lengths of 
multi-loop X-class flares derived using the R97 approach can be comparable to the observed 
loop lengths, considering that the application of the single-loop approach to multi-loop flares 



1 Images of complex solar flaring loop arcades spatially resolved by TRACE ( Handv et al at Extreme 

Ultraviolet (EUV) wavelengths and their associated X -ray light curves obtained by the Geostationary Oper- 
ational Environmental Satellites ( GOES, IGarcialll994 and references therein) can be found and inspected at 
the TRA CE flare catalog web site |http : / /hea-www . harvard . edu/trace/f lare_catalog/index . html (see 
also Figure IT2|). 
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is a priori incorrect. 

The consistency may be by chance. For instance, for all 5 flares, the inferred slope ( 
in the logT - log V EM diagram is found to be close to ~ 0.4. Let us ignore for a moment 
the fact that the physical meaning of the slope £ for multi-loop flare arcades is generally 
different from that of a single loop (i.e. different plasmas in different loops heated and 
cooled sequentially versus a single plasma in a single loop heated and cooled). Within the 
framework of the single- loop model of R97, ( ~ 0.4 is close to the lowest values allowed, and 
the equations of R97 produce smaller loop lengths at smaller slopes. In fact, in the regime 
of low ( (C < 0.7), the loop length is steeply decreasing with decreasing (, by a factor of 10 
when ( changes from 0.7 to 0.4. 

Contrary to the "by chance" explanation above, loop length numbers can be comparable 
due to the presence of a single loop (or localized multiple loops ignited simultaneously and 
undergoing similar heating and cooling processes) that dominates the flare X-ray emission. 
For instance, an arcade-like structure with a single primary loop dominating the rise and 
early decay phases of a flar e was proposed (b ased on a detailed modeling) for the complex 



flare on Proxima Centauri ( jReale et al.l 120041 ) . At least for the best studied of the 5 solar 



X-class flares considered here, the Bastille Day flare, we do no find compelling reasons for 
the presence of such a dominant X-ray emitting structure (Appendix |A]) . 

Clearly, even if the loop lengths inferred via the single-loop method happen to be com- 
parable to the observed values for such complex flare arcade events, the method would not 
predict a correct loop geometry and would likely underestimate the emitting volume by at 
least a factor comparable to the number of instantaneously heated loops in the arcade, Ni oop , 
and would overestimate the average electron density of emitting plasma by ~ \jNi oop . For 
instance, during the Bastille Day flare arcade, Ni nnp ~ 20 out of the ~ 100 observed loops 



fired near the flare peak time (lAschwanden Alexander! l200ll ) 



For all 5 solar X-class flares considered here, the inferred slope ( in the log T - log \JEM 
diagram is found to be close to the val ue of ~ 0.4. Observati onaUy, such a value is among the 



lowest values seen in solar flares (e.g. ISylwester et al.lll993l ). Such a shallow slope indicates 



continued heating during the decay phase of the X-ray emission. Continued heating is also 
what one would expect for multi-loop two-ribbon flares where the decay phase could be 
entirely driven by the heating released through sequential reconnection in individual loops 



flKopp fc Polettd 11984 1. 

The slope ( > 1 in the log T - log y/EM diagram for the DQ Tau flare is significantly 
larger than the ( ~ 0.4 found for complex powerful solar flares (Appendix |AJ , suggesting 
heating behaviour different from the 5 solar X-class flares. The possibility that the DQ Tau 
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flare could involve an arcade-like structure with a single loop dominating the rise and earl y 
decay phases, as was proposed for the complex flare on Proxima Centauri (IReale et al.ll2004j ). 
can not be excluded. On the other hand, within th e framework of a model involving explosive 
chromospheric evaporation into a single loop (e.g. iReald 120071 ). an evaporation time of only 
3 — 5 ks is needed for a plasma with temperature T Q b s ~ 60 — 90 MK ( §3.3[) to propagate with 
an isothermal sound sp eed v s [cm s _1 ]= \J (5/3kbT b s )/ {^m p ) = 1.5 x 10 4 -^/(r o 6 S [K]) (e.g. 
Aschwanden et al.l 120001 ) and to fill a DQ Tau flaring loop of a semi-length L = 4 — 5 R*. 
The observed rise timescale of the DQ Tau light curve, r rise = 26.3 ks, is much longer than 
3 — 5 ks and thus might indicate a complex history of heating during the early phase of the 
flare, e.g., due to multiple loops. 

It is also worth mentioning that the time delay between temperature and emission 
measure peaks observed in the DQ Tau flare ( §3.3p can not be used to argue in favor of a 



single loop event. This effect is commonly seen in multi-loop flares and can be explained by 
invoking the principle of linear superposition (Appendix |XJ). Likewise, the Neupert-like effect 
observed in the DQ Tau flare (§!]) can not be used to argue in favor of a single loop event. If 
the Neupert effect occurred in each loop of a multi-loop system, then an observation that did 
not spatially resolve the loops would show the Neupert effect. Low-resolution obs ervations 
of multi-loop solar flares commonly exhibit the Neupert effect (jVeronig et all 120021 ) . 



After all, in view of the fact that the observed timescales, temperatures, and X-ray 
luminosities of PMS stars (DQ Tau and ONC stars) are much higher than those of typical 
solar flares (§5]), one could argue that the solar analogy might not be directly applicable. 
PMS and magnetically active older stars possess stronger, than solar, surface and global 
magnetic fields and larger volumes for magnetic fields to interact, so the ext reme flaring 
behaviour, beyond solar analogy, is expected (e.g. §2.2 in lBenz fc Giidelll2010l . and §6.21 in 
this work). But then the single- loop approach developed on solar analogy might not be 
applicable either. 

We conclude that we can neither prove the presence of loop arcades analogous to the 
solar cases, nor refute the presence of single flaring loops longer than any seen on the Sun. 
The results from the single-loop X-ray modeling presented in §3.41 and §4.11 should not be 
treated as definitive and should be considered with caution. 



6.2. Origin of the DQ Tau Flare 



The origin of the big flares from the COUP sample themselves is unclear. For some 
of the 32 most powerful COUP flares whose inferred coronal structures reach several stellar 
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radii, iFavata et al.l (120051 ) suggest magnetic loops linking the stellar photosphere with the 
inner rim of the circumstellar disk. For the extended sample of > 200 flares G08ab find 
that large (a few to several stellar radii in length based on the single-loop approach of 
R97) coronal structures are present in both disk-bearing and disk-free stars; however G08ab 
also find a subclass of super-hot flares with peak plasma temperatures exceeding 100 MK 
that are preferentially present in highly accreting systems. G08ab further propose that the 
majority of big flares from the COUP sample can be viewed as enhanced analogs of the 
rare solar long-decay events (LDEs; see §7.3 in G08b and references therein). LDEs are 
eruptive solar flare events that produce X-ray emitting arches and streamers with altitudes 
reaching up to several hundred thousand kilometers (L ~ 0.5 — 1 Rq). In panels (b, d, e) of 
Figure [101 representative solar LDEs (gray diamonds; compiled in G08a) are shown to have 
systematically higher flare durations and coronal lengths than the more typical solar flares 
(solid and dotted gray loci). We wish to emphasise that the analogies of the big flares from 
the COUP sample to solar LDEs ( "growing" systems of giant multiple loops) are not based 
on information about the specific detailed geometry of COUP flares (the geometry is really 
un- known to us and is only simplistically modeled as a large single loop), but are instead 
based on the simple fact that the observed flare durations and model-inferred characteristic 
loop scales for COUP flares are the largest ever reported from PMS stars, provided that these 
model-inferred loop scales are close to the truth. The most widely accepted model for the 
origin of LDEs is that the impulsive flare near the solar surface (L < 10 -2 R ) blows open 
the overlying large-scale magnetic field with subsequent reconnection of large-scale magnetic 
lines through a vertical current sheet. The large-scale magnetic field of PMS stars is likely 
far stronger than in the Sun, and can sustain giant X-ray arches and streamers with sizes 
L ~ 1 — 10 R*. For the DQ Tau flare a different mechanism from that of a solar LDE-like 
reconnection can be considered. 

There are at least three independent supporting lines of evidence suggesting that the 
DQ Tau flare could be produced through a process of colliding magnetospheres. 

1. Although the DQ Tau flare duration, morphology, and plasma temperature are typical 
of those of big flares from the COUP sample (§5}, the probability of observing a 
big COUP-like flare is small. Within the COUP observation window (~ 1200 ks) 
on average 3 big flares per bright PMS star were detected. This number is derived 
from the analysis of the 161 brightest X-ray PMS stars (G08a); the frequency of big 
flares for the remaining 1200 fainter COUP PMS stars is even lower. The DQ Tau 
Chandra observation was designed for an X-ray flare to appear close to the orbital 
phase of $ = 0.98, wh ich corresponds to the peak of the April 2008 large sub-mm 



flare (ISalter et al.l 120081 ). The detected Chandra flare indeed peaks within 15 ks of 



this orbital phase. The probability of detecting a big COUP-like flare within 15 ks of 
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a pre-specified point in time is small (Prob < 15 x 2/(1200/3) = 0.075). 



Bower et al. 


2003; 


— „„„ ^ 
Salter et al. 


2010) 



in PMS stars has been rarely reported (e.g. 
DQ Tau X-rays are accompanied by a relatively unique non-thermal mm activity. The 
X-rays and mm are likely related through a Neupert-like effect O- The observed 
re-appearance of this unique mm activity near periastron is proposed to be associa ted 
with processes caused by the interacting magnetospheres (ISalter et al.ll2008l 120101 ). 



3. It is natural to expect large-scale magnetic structures to be involved in a process of 
colliding magnetospheres. The loop flaring sizes inferred from the X-ray flare analysis 
(4 — 5 R*) are comparable to half the separation of the DQ Tau components near 
periastron ( §3.4.1[) . This third item of evidence is completely model-dependent and 
should be considered with caution ( §6.ip . 



6.3. Energetics and Loop Geometry 



Based on the concept of interacting magnetospheres, we can further comment on ener- 
getics and possible flare loop geometry. 



We can show that a crude estimate of the power expected to be dissipated by magnetic 
reconnection from colliding magnetospheres in DQ Tau is in agreement with the derived peak 
flare X-ray luminosity (£x,o.5_8 = (4 — 5) x 10 30 ergs s -1 ; see Table |2]) and the modeled rate 
of the kinetic energy injected into the chromosphere by non-thermal electrons (a x Fn(ti) = 
(0.6 - 0.8) x 10 32 ergs s _1 ; gj). 

For a large scale dipolar topology and a magnetic field B(L) ~ B ph j{Lj R+ + l) 3 with 
photospheric field strengths in the range 1 — 6 kG consistent with measurements of Zeeman 



broadening and circular polarization of photospheric lines in PMS stars (e.g. iJohns-Krull 
20071 : iDonati et al.ll2008l ). the field strength at distances of 4 — 5 R* from the stellar surface 
is expected to be B(4 — 5R*) ~ 5 — 50 G. Within the framework of the simple magnetic 
reconnection model shown in Figure [TTJ the initially dipolar fields of average st rength B = 5— 
50 G at 4 — 5 R* from the stellar surface colliding at a speed of v ~ 100 km/s (IMathieu et al. 



19971 ) reconnect with the rate of energy release E m /rR ~ 10 30 — 10 32 ergs s _1 , which is 
in agreement with the numbers given above. It is interesting to note that the analogous 
procedure applied to cases of star-planet magnetic interaction have difficulty explaining an 
excess X-ray emission associ ated with st ellar chromospheric hot spots rotating synchronously 
with close-in giant planets (ILanzal 120091 ) . 
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The model of interacting magnetospheres in RS CVn-type binaries by lUchida fc Sakurai 



(119851 ) predicts the existence of X-ray emitti ng loop struc tures connecting the stars. The 
model of star-planet magnetic interaction by lLanzal ( 120091 ) expects complex topologies, in- 
cluding loops connecting the planet with the stellar surface. While our X-ray modeling 
method and many related results are limited to the scenario of a single X-ray emitting loop 
on a single star, other scenarios for the DQ Tau loop geometry are possible: multiple X-ray 
emitting loops on a single star ( §6 ■ 1 j) . two or more loops appearing on both stars simultane- 
ously, or a single or multiple loops connecting the two stars. For our modeled scenario, where 
the loop's footprints are anchored on the surface of a single star, potential destruction of the 
loop by centrifugal forces is not a serious issue. The derived length for the X-ray emitting 
structure of 4 — 5 R* is only comparable to and does not exceed the Keplerian corotation 
radius R cor =4.7 rM 



Conclusions 



The young high-eccentricity binary DQ Tau exhibits po werful recurring mm flari ng at- 
tributed to collisions between the two stellar magnetospheres (jSalter et al.ll2008l |2010| ). The 
separation of only ~ 8 R* between the binary components at periastron implies magneto- 
spheric interactions at characteristic scales and magnetic strengths comparable to those of 
relatively rare big X-ray flares from single COUP PMS stars (G08ab). To search for X- 
ray emission arising from processes associated with colliding magnetospheres, we performed 
simultaneous Chandra X-ray and IRAM/CARMA/SMA mm observations of an orbital seg- 
ment of DQ Tau near the peaks of the previously detected mm flares. We report here a 
Chandra detection of a long X-ray flare accompani ed by mm ac t ivity close to periastron. 
The mm observations are discussed in more detail in ISalter et al.l (120101 ) . 



We start by addressing light photon pile-up in the Chandra data ( §2.ip . Further com- 
parison with the constant and soft X-ray emission of DQ Tau at an orbital phase away 
from periastron detected by XMM-Newton shows similar spectral characteristics between 
the time-integrated XMM-Newton spectrum and the low temperature plasma component of 
the time- integrated Chandra spectrum ( §3.21) . Derived from the single- loop flare model of 
R97, the half-length of the X-ray emitting structure of L = 4 — 5 R* is comparable to half 
the binary component separation at periastron ( §3.4.11) . Inferred from the cooling time-scale 



15 The Keplerian corotation radius where the centrifugal force balances the gravitational force is R cor — 
(G x M x P 2 /4 x 7T 2 ) 1 / 3 . For the DQ Tau binary components with similar masses of M = 1.6 M Q and 
rotational periods of P — 3 days, R cor = 4.7 R*. 



-25 - 



analysis, the cross-sectional radius of the loop is r = 0.03 — 0.06L ( §3.4.2(1 . In view of a pos- 
sible solar analogy, the concept of the single-loop approach and its application to DQ Tau 
flare might be questionable ( §6.111 . Thus, both results on the loop length and cross-sectional 
radius of the loop should be considered with caution. The coincidence of the X-ray flare peak 
with the end of one of the mm flares, and the resemblance of the thermal energy profiles 
derived independently from the microwave and X-ray single-loop modeling, give an indica- 
tion for a Neupert-like effect (jjlj). Comparison with the large flares from the COUP sample 
shows that the DQ Tau flare duration, morphology, and plasma temperature are typical of 
COUP flares (35}. Three lines of evidence (the low probability of detecting a COUP-like 
flare, the Neupert effect apparent in the relative timing of the X-ray and mm flares, and the 
similarity between the binary separation and the inferred loop size) support the notion that 
the DQ Tau flare could be produced through a process of colliding magnetospheres ( §6.2p . 
The third evidence is fully model-dependent and should be considered with caution ( §6.ip . 
The inferred flare peak X-ray luminosity and the modeled rate of the kinetic energy injected 
into the chromosphere by non-thermal electrons are in agreement with an estimate for the 
power dissipated by magnetic reconnection within the framework of a simple model of inter- 
acting magnetospheres ( §6.311 . Our flare analysis is limited to the scenario of a single X-ray 
emitting loop on a single star; other scenarios for the DQ Tau loop geometry are possible: 
multiple X-ray emitting loops on a single star ( §6. 1(1 . two or more loops appearing on both 
stars simultaneously, or a single or multiple loops connecting the two stars ( §6.3p . 

Future coordinated multi-wavelength observation campaigns of DQ Tau, spanning a 
larger range of orbital phase and more wavebands than in this study, are highly desirable. 
These will allow further investigation of the physical processes during an interaction of 
magnetospheres. Due to the predictable and bi-weekly recurrence of flare activity near peri- 
astron, DQ Tau offers an excellent opportunity to persistently study the poorly understood 
connection between particle acceleration and coronal heating in stellar flares. 
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A. Application Of The Single Loop Model To Solar X-class Flares 



R97 applied their single loop approach to solar X-ray flares observed with Yohkoh/SXT 
and showed that the method provided loop scaling sizes comparable to those measured in 
SXT images. In this section we further check on the validity of the method by applying 
it to solar X-class flares clearly associated with multiple loop arcade geometries spatially 
resolved by TRACE. To better mimic stellar X-ray observations, which cannot resolve flare 
structures, the analysis of these solar X-class flares is carried out on disk-integrated X-ray 
data obtained by GOES. 



From the TRACE solar flare catalog 16 ! we choose only limb flares, where the heights of 
coronal structures derived from high resolution TRACE images at EUV are most reliable. 
Out of the 75 X-class flares listed in the catalog, ~ 10 are limb flares, and of those at 
least the following 4 have been studied in detail: X3.7 and X2.5 flares o n November 22, 
1998 terren||2QQ0); X1.5 flare on April 21. 2002 (IcTallagher et al1 l2002h: and X3.1 flare 
on August 24, 2002 Li fc GarJl2005l : Ifleznikova et allboOflf l. To these 4 flares we add the 
well studied Bastille Day (14 July 2000) near-disk center flare (lAschwanden fc Alexander 
200ll ). All these five flare events are ass ociated with arcades o f multiple loops ( Figure [T2l . 
At le ast for the two events on 21-Apr-02 jGallagher et al.ll2002h and 24-Aug-2002 \u fc Gan 
20051 ). "rising" flaring loop systems are reported and are believed to be associated with the 
formation of new loops as magnetic reconnection progresses upward. The heights of the 
EUV loops inferred from TRACE images, H, are given in Column 9 of Table |3J We assume 
that the X-ray emitting loops for the five events have heights comparable to or higher than 
the observed EUV loops; this assumption is generally expected d ue to effects o f loop cooling 
and shrinkage, and is found for both 22-Nov-98 flares (Figure 3 in I Warren! 120001 ) . As a sanity 
check of our GOES X-ray flare analysis, presented below, the 3 brightest M-class flares from 
the original testbed sample of R97 are also analysed. 



The GOES X-ray data have been analysed using the GOES software tools within an 



fottp : //h ea-www . harvard . edu/ trace/f lare_catalog/ index . html 
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IDL based data analysis software package for Solar Physics, SolarSoft (IFreeland fc Handy 
1998ll 17 l . Background subtraction for these powerful solar flares is ignored. Plasma tem- 



perature and emission measure pr ofiles are derived by applying the filter-ratio method 
(goes-chiantLtem.pro program from IWhite et al.l 120051 ) to the observed GOES light curves 
in the (0.5 - 4 A) and (1 - 8 A) bands. 

Light curve timescales, observed peak plasma temperature and emission measure, and 
slope £ of the trajectory in the log T - log V EM diagram are shown in Figure d2] and are given 
in Columns 2 — 6 of Table [21 Estimates on loop half-length from the rise, L rise , and decay, 
Ldecay, phases of light curves, using the single loop approach, are also tabulated (Columns 7— 
8) and compared to the heights of the corresponding EUV loops, H (Column 9). Loop 
half-length from the rise timescale of the light-curve, assuming an explosive chromospheric 
evaporation, is given as L rise [cm] = r rise [s] x 1.5 x 10 4 A/T obspk [K] ( §6.ip . Loop half-length 
from the decay phase of the light-curve is calculated using the R97 method ( §3.4. ip with 
the instrument-dependent para meters in th e calibration formulas for F(() and T p k taken 
for GOES 9 from Table A.l of iReald ( 120071 ). Formal statis tical errors on der ived GOES 
temperature and emission measure are expected to be < 4% (1 White et al.l 120051 1 . leading to 
statistical errors on the slope ( of < 3% and errors on the loop half-length derived from the 
decay phase of the light curve L decay of < 30%. 

For the three brightest M-class flares from the original sample of R97, comparison of 
their flare properties derived using Yohkoh/SXT data (R97) with those inferred from our 
GOES analysis shows that GOES peak plasma temperatures are systematically higher by 
30%, slope ( differences are within 30%, and loop half-length L decay differences are within 
50% (with a formal statistical precision of < 30% on individual flares). Effects of neglecting 
background emission in GOES data are expected to be more pronounced for M-class rather 
than for X-class flares. 

The result of our analysis is that for the five X-class flares associated with arcades of 
multiple loops (Figure [12]), loop half-lengths Ldecay (or more precisely loop heights, 0.64 x 
Ldecay, assuming semi-circular loops stand vertically on the solar surface) derived from disk- 
integrated GOES data using the single loop approach of R97 are shorter or comparable to the 
heights of the EUV loops, H, measured from a direct inspection of high- resolution TRACE 
images (Table [3]). It is important to note that the application of the single-loop approach to 
multi-loop flares is a priori incorrect, and the similarity in loop length values can occur by 



17 Descriptions and codes for SolarSoft, and GOES IDL User Guide can be found at 
|http : // www.lmsal . com/solarsof t/ and http://hesperia.gsfc .nasa.gov/$\sim$kim/goes_sof tware/goes . html| 
respectively. 
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chance ( §6.11) or because of a special physical phenomenon taking place, namely the presence 
of a dominant X-ray emitting structure. The latter does not seem to apply to the Bastille 
Day flare (see below). Meanwhile, the loop half-lengths derived from the rise phases of the 
GOES light curves with the assumption of chromospheric evaporation in a single loop, L rise , 
(or their semi-circular heights, 0.64 x L r i se ) are systematically larger by a factor of 4 — 7 
compared to the measured heights of the EUV loops. It is also important to note two other 
aspects of our analyses. First, for all five flares the inferred slope of the trajectory in the 
logT - log y/EM dia gram is ( ~ 0.4, indicating that heating is the dominant process in 
the decay phase of a flare ( §6.11) . Second, the effect of time delay between temperature and 
emission measure peaks is clearly observed in all 5 multi-loop X-class flares (Figure [121) . This 
confirms the general expectation that the effect (described in §3.31) is applicable to a multi- 
loop flare case as well as to a single loop. In a multi-loop case the effect can be explained 
by the principle of linea r superposition of flares from ind ividual loops to produce a flare for 



a complete event (e.g., lAschwanden fc Alexander! l200ll ). Below we apply the principle of 



superposition to the Bastille Day flare. 

If individual loops in a flaring loop system are ignited nearly simultaneously and are 
subject to a similar evolution, which is rather a quite rare condition, (or a single loop domi- 
nates the flare X-ray emission), then for the complete event, according to the superposition 
principle, the flare temperature profile (including its peak plasma temperature), the flare 
light curve shape (specifically the e- folding decay timescale), and the slope ( on the logT - 
log VEM dia gram should be similar to those of individual loops. In such case, for the 
complete event the equation of R97 would yield loop length similar to those of individual 
loops (or a single dominant loop). Is this the case for the 5 solar X-class flares studied here? 
At least for the best studied case, the Bastille Day flare, we are unable to collect strong 
evidences in support of such a picture. 

The plasma temperature profile of the Bastille Day event (Figure [T2"j) has a major 
peak (near 10:18:20 UT) and two bumps seen along the rise (near 10:10:00 UT) and the 
decay (near 10:28:20 UT) phases of the profile, respectively. The classic time delay from 
temperature to emission measure is seen at the major peak (AT ~ 340 seconds) and 
decay-phase bump (AT ~ 200 seconds). According to the TRACE images (Figure 7 in 
Aschwanden fc Alexander! l200ll ) and the movi the rise-phase temperature bump is likely 



related to the appearance of several EUV loops at the westernmost edge of the arcade (near 
10:10: 00 UT), the major temperature /e mission measure peaks are associated with ~ 20 EUV 



loops (lAschwanden &: Alexander! l200ll ) that start appearing at about 10:14:00 UT (and last 



3 The TRACE movie of the Bastille Day flare can be found e.g., at 



http : //www. suntrek. org/gallery/gallery-magnetic-sun. shtml 
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up to 11:00:00 UT), and the later temperature/emission measure bumps are likely related 
to numerous loops at the eastern part of the arcade that produce weaker flares appearing 
around 10:38:00 UT. 

Using the principle of superposition, we simulate profiles of the complete event as 
EM tot {t) = J2 EM i(t) and Ttotif) = (Tiit) x EMi{t))/Y,EMi{t) where the summation is 
over the 20 EUV loops associated with the major peak. Here we assume that the Ti(t) and 
EMi(t) profiles for individual loop events have shapes similar to the profiles of the observed 
complete event (Figure [T2"|) . This simplistic simulation shows that if the loops fire nearly 
simultaneously (within the time period less than the observed T — EM peak delay of 340 sec- 
onds, or even more precisely less than 2 minute period), then the temperature profile (both 
the shape and the peak), the shape of the emission measure profile, as well as the T — EM 
peak delay of the complete event look similar to those of individual loop events. Otherwise, 
for a complete event the simulation predicts broader shapes of the temperature and emission 
measure profiles and longer T — EM peak delay than those of individual loop events, e.g., the 
delay becomes twice longer, if the loops fire sequentially (and evenly spaced in time) within 
the 10 minute time period. Thus, the ignition of 20 loops (that are subject to a similar evolu- 
tion) within the 2 minute time period might mimic a single dominant loop. Notice, however, 
that in our analysis we make an assumption that the profiles for individual loop events have 
shapes similar to the observed profiles of the complet e event. Many other solutions mu st ex- 
ist. For example, using the principle of superposition lAschwanden &: Alexander! ( 1200 ll ) show 
that individual loop events not firing near-simultaneously, with temperature and emission 
measure profiles much shorter than for the complete event, can produce the observed profiles 
of the complete event. 

One might also argue that EUV images might not provide a true picture of an X-ray 
flare event. The X-ray band Yohko h/SXT image of the Bastille D ay flare corresponding 
to the major EUV peak (Figure 3 in lAschwanden fc Alexander! 120011 ) shows a single bright 
feature (spanning 20 x 40 square arcseconds) within the more extended weaker emission, 
which may be eyidence that a few localized EUV-like loops^l dominate the X-ray emis- 
sion (IReale et al.ll2004f). Howe ver, since only a single instant X-ray image (Figure 3 in 
Aschwanden fc Alexander! 1200 ll ) is available, it is unclear how the X-ray emission evolves 
spatially during the early phase of the flare. 

Assuming the presence of a dominant X-ray emitting structure (a single loop or multiple 
similar loops firing simultaneosly) during the Bastille Day flare we can apply the R97 ap- 



19 In the EUV data (Figure 8 in lAschwanden fc Alexanderll200ll ) we estimate that a few to several EUV 
loops can fit within the apparent footprint of the bright X-ray feature. 
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proach strictly to the initial decay segment of the X-ray light curve (IReale et al.ll2004[ ). The 
slope of the initial decay (prior to the bump at the late decay phase) on the log T - log V EM 
diagram is ( ~ 2, much steeper than the ( = 0.41 obtained for the entire decay phase of the 
light curve, leading to an unrealistically large loop size of 200,000 km, which is a factor of 
10 larger than the loop lengths measured directly from the EUV images. Considering the 
above, it seems to us that there is no clear evidence in support of a dominant X-ray emitting 
structure in the Bastille Day flare. 
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Table 1. X-ray Models for the Low Temperature Components 



Parameters 


Model 1 


Model 2 


(i) 


(2) 


(3) 


N H (10 21 cm" 2 ) 


1.3 


3.2 


fcTi (keV) 


0.7 


0.3 


EM X (10 52 cm -3 ) 


4.1 


19.1 


kT 2 (keV) 


1.9 


0.9 


EM 2 (10 52 cm -3 ) 


6.6 


8.3 


L x (10 30 org s- 1 ) 


1.0 


1.7 



Table 2. Time- Resolved Spectroscopy 



Time Segment Using Model 1 Using Model 2 
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Note. — Columns 1-3: Segment's sequential number, central time, and duration. Columns 4-7: Using Model 1 for the fixed low temperature plasma component 
emission, the flare properties inferred from the spectral fits arc: the flare plasma temperature and its lcr errors, emission measure and its lcr errors, absorption- 
"corrcctcd" luminosities in the (0.5 — 8) and (0.01 — 50) kcV bands. Columns 8-11: Same as in Columns 4-7, but using Model 2 for the low temperature plasma 
component emission. 
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Table 3. Inferred Properties of The Testbed Solar Flares 



Flare 


T rise 




Slope 




EM k 




Ldecay 


H 


Reference 




(s) 


(s) 




(MK) 


(10 49 cm" 3 ) 


(1000 km) 


(1000 km) 


(1000 km) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


22-Nov-98 06:30 X3.7 


77 


515 


0.38 


23.0 


19.6 


56 


1 


£ 10 


1 


22-Nov-98 16:10 X2.5 


119 


655 


0.48 


22.4 


13.2 


84 


11 


> 10 


1 


14-Jul-00 10:00 X5.7 


213 


1422 


0.41 


19.3 


29.8 


140 


8 


17.5 


2 


21-Apr-02 00:05 XI. 5 


761 


3889 


0.41 


15.5 


9.4 


453 


18 


30 - 50 


3 


24-Aug-02 00:30 X3.1 


215 


1680 


0.42 


23.0 


17.8 


154 


14 


> 15 


4, 5 


02-Nov-91 06:44 M9.1 


195 


737 


0.55 


20.5 


3.2 


134 


18 




6 


06-Feb-92 03:16 M7.7 


281 


955 


0.93 


19.8 


2.5 


189 


64 




6 


15-Fob-92 21:29 M5.5 


277 


712 


0.83 


17..-) 


2.1 


174 


38 




6 


Note. — Columns 1: 


Flare date 


start 


time, and 


GOES class. 


Columns 2-3: 


Flare light-curve rise and det 


:ay timcscalcs. 


Column 4: 


Inferred slope of the trajectory in 


the log 


T - log V 


EM diagram. 


Columns 5-6: 


Flare peak 


bserved plasma 


temperature and emission 



CO 



measure. Column 7: Loop half-length from the rise timcscalc of the light-curve, assuming an explosive chromosphcric evaporation in a single 
loop, L r j se [cm] — r r i sc [s] X 1.5 X 10^ ^/^obs.pk P^] • Column 8: Loop half-length from the decay phase of the light-curve using the single loop 
approach of R97. Column 9: Apparent loop height from TRACE images. Column 10: Reference for columns (1) and (9). REFERENCES: 
(D lWarrer] j 20001) ■ (2) lAschwanden &: Alexanderl teOOH ), (3) iGallagher et al.| j2002D . (4) |Li & G"H i2005l) . (5) iReznikova et all i2009l) . (6) 
R97. 
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Time [ks] 

Fig. 1. — The January 2010 59 ks Chandra observation of DQ Tau at the orbital phase $ = 0.95 — 0.99 close 
to periastron passage $ = 0.0, 1.0. (a) X-ray events detected in a 1.4' x 1.4' field around DQ Tau. Source 
extraction aperture enclosing 99% of the local PSF is marked by the solid "flower" contour. Background 
extraction region (dashed annulus) excludes > 99% of the PSF power. Direction of the faint CCD readout 
streak emanating from DQ Tau is indicated by the arrowheads, (b) Observed light curves of DQ Tau 
binned by 1 hr in the full 0.5-8 keV (open circles), soft 0.5 — 2 keV (triangles), and hard 2 — 8 keV 
(squares) energy bands. Light curve in the full band, corrected for mild photon pile-up, (filled circles; £)2.2[) 
with the best-fit exponential rise and decay (solid lines) models (e-folding timescales of T r i se = 26.3 ks 
and Tdecay = 40.9 ks, respectively), (c) Evolution of the adaptively smoothed median X-ray event energy, 
observed (open circles) and corrected for pile- up (filled circles). The smoothing is performed using a sliding 
boxcar kernel enclosing 1000 counts. For the six independent time segments (segments ## 1, 11, 22, 27, 36, 
45 from Table [5]) normalized median absolute deviation (MAD) errors on median energy (see Appendix B 
in G08a) are indicated by the vertical bars. 
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Fig. 2.— The February 2007 12.6 ks XMM-Newton observation of DQ Tau at the orbital 
phase $ ~ 0.65 - 0.7. (a) Merged EPIC MOS+PN (0.2 - 10) keV image of the ~ 15' x 15' 
region around DQ Tau. The extraction aperture (small solid circle) encloses 90% of the local 
PSF power. MOS and PN background extraction regions are marked by the large and small 
dashed circles, respectively, (b) Background subtracted EPIC-PN light curve of DQ Tau in 
the (0.2 — 10) keV energy range with a bin size of 100 s. 
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Fig. 3. — Comparison of the time-integrated Chandr a- ACIS and XMM-Newton-EPIC spec- 
tra, and the ambiguity in the spectral fitting. Panels (a) and (c) show spectral fits to the 
Chandra and XMM- Newton data, respectively, with their Nh, kT%, kT2, EM2/EM1 param- 
eters fixed at values of the low temperature plasma component Model 1 (Tabled]). Chandra 
modeling requires an additional time-averaged flare component with kT 3 = 4.3 keV. Pan- 
els (b) and (d) show the Chandra and XMM-Newton fits, respectively, with Nh, kTi, kT-i, 
EM2I 'EM\ fixed at values of the low temperature plasma component Model 2 (Table [I]). 
Chandra modeling requires an additional time- averaged hot temperature (flare) component 
with kT 3 = 3.5 keV. XMM-Newton modeling requires an additional moderately hot compo- 
nent kT 3 = 1.9 keV. 
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Fig. 4. — Fits of the pile-up corrected Chandra spectra for the six independent time segments 
(segments 1, 11, 22, 27, 36, 45 from Table [2]) exemplify our time-resolved spectral 
modelling. For each of the segments the 2-T soft plasma component model is frozen (two 
leftmost dotted lines) while the plasma temperature and emission measure of the 1-T hot 
(flaring) component (rightmost dotted line) are varied to obtain the best-fit. Modeling is 
shown for the case of the low temperature plasma component Model 1 (Tabled]). 
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Fig. 5. — Evolution of the inferred plasma temperature (a) and emission measure (b) of the 
Chandra DQ Tau flare using low temperature plasma component Model 1 (see Table [2]). (c) 
Evolution in the log T - log VEM plane with the best-fitting flare decay as a dashed line 
of slope £■ Black (and gray) points are shown for the cases of the pile-up corrected (and 
uncorrected) data. Temperature and emission measure points are derived for overlapping 
time segments ( §3.3[) . thus their errors are not independent. 
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Fig. 6. — Evolution of the inferred Chandra DQ Tau flare properties using low temperature 
plasma component Model 2. See Figure [5] for details. 
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Fig. 7. — Comparison of the derived flare cooling timescales with the observed light curve 
decay timescale for a testbed X-ray flare, specifically the superflare from COUP source 
# 1343. (a) Evolution of the inferred plasma temperature (solid line), emission measure 
(dotted lines), and absorption-corrected X-ray luminosity (dashed line). The flare decay 
phase of interest is indicated by the thick line, (b)-(f) Evolution of the derived e- folding 
cooling time scales: thermal conduction cooling (r con , thin dashed line), radiation cooling 
(r ra( i, thin solid gray line), combined conduction and radiation cooling r t h iXl r th = ^/Tcon + 
t/ T rad, thin solid black line). The e-folding decay timescale of the observed Chandra light 
curve (Tdecay) is indicated by the thick solid line. Cooling time scales are derived assuming 
a flaring loop length of L = 4.8 R* and five different values for the ratio of the loop cross- 
sectional radius to the loop length (3. 




Fig. 8. — Comparison of the derived flare cooling timescales with the observed light curve decay 
timescale for the Chandra flare toward DQ Tau using low temperature plasma component Model 1. 
(a) Evolution of the inferred plasma temperature (solid line) , emission measure (dotted lines) , and 
absorption-corrected X-ray luminosity (dashed line). The flare decay phase of interest is indicated 
by the thick line, (b)-(i) Evolution of the derived e- folding cooling time scales: thermal conduction 
cooling (r con , thin dashed line), radiation cooling (r ra ^, thin solid gray line), combined conduction 
and radiation cooling (1/t#, = l/r con + 1/T ra d, thin solid black line). The e- folding decay 
timescale of the observed Chandra light curve is indicated by the thick solid line. Panels (b)-(e) 
present cooling timescales derived assuming the loop length of L = 5 R* for four different values 
of (3. Panel (c) here shows the timescale pattern reminiscent of that seen in COUP # 1343 flare 
(Figure [7b). Panels (f)-(i) exemplify evolution of cooling timescales for loop lengths both less 
than and higher than an L = 5 R*, with the values of the loop thickness parameter f3 selected to 
emphasize cases that most closely follow the cooling timescale evolution pattern seen in panel (c). 
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Fig. 9. — (a) Comparison between the IRAM mm (small gray circles) and Chandra X-ray (big 
black circles) light curves for DQ Tau. The mm flare decay e-folding timescale is tr = 15.9 ks 
(dashed line), (b) Indication of a Neupert-like effect. Plasma thermal energy (thick solid line) 
inferred from our X-ray modeling (using Model 1) compared to the convolved IRAM light curve 
using formula (12) with parameters C\ = 0.5 x 10 35 ergs and a x Fn(ti) = 0.8 x 10 32 ergs s _1 (thick 
dashed line). In the X-ray modeling, the energy points are derived from overlapping time segments 
(£ )3.3|) ; thus the energy errors are not independent. 
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Fig. 10. — Comparison of the DQ Tau flare properties with those of the COUP PMS stars 
(panels (a)-(e)) and flares on older stars (panels (b), (d), (e)). (a) Flare decay vs. rise 
timescales. (b) Flare duration vs. peak flare temperature, (c) Time-integrated flare energy 
(determined from the flare peak X-ray luminosity times the flare decay timescale) vs. peak 
flare temperature, (d) Peak flare emission measure vs. peak flare temperature, (e) Flare 
loop length vs. peak flare temperature. Large flares from the COUP sa mple ar e show n as 
black circles. In panels (b), (d), and (e), the stellar flares compiled by iGiidell (120041 ) are 
marked by gray boxes with their linear regression fit and la range (in panel (d)) shown as 
gray dashed-dotted and dashed lines, respectively. Regression fit and la range for solar- 



stellar (panel (b)) or just solar (panels (d) and (e)) flares compiled by lAschwanden et al. 



( 120081 ) are marked by gray solid and dotted lines with representative solar LDEs (compiled 
by G08a) discussed in the text shown as gray diamonds. The DQ Tau flare is marked by the 
thick black square with bars representing systematic differences between the flare properties 
derived using modelling with the low temperature plasma component Model 1 and Model 2. 
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Fig. 11. — Illustration of the process of interacting magnetospheres in DQ Tau (in 2D 
projection plane). Just before periastron, with a component separation of ~ 10 R*, the 
magnetospheres collide with velocity v ~ 100 km s _1 . Large-scale magnetic field lines of 
both components at > 4 R* above the stellar surface violently interact (lines above 5 R* are 
ignored here). On the plane perpendicular to the velocity vector, the interacting magnetic 
lines occupy an area A ~ 6 x B% (hatched; assuming dipolar topology). The total amount 
of magnetic energy within the interacting volume of width / for both stars is E m ~ 2 x 
(B 2 /8ir) x I x A and this is released over a timescale tr~1/v. 
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Fig. 12. — Modeling of five solar X-class flares, each represented by a single row of the 
figure. TRACE images (171 A or 195 A filter) are shown in the first column; fields of view 
are 100" x 100" (panel a), 100" x 100" (panel e), 280" x 128" (panel i), 210" x 210" (panel 
m), and 110" x 110" (panel q). GOES light-curves (1 — 8 Aband, gray) with least-squares 
fits of the rise and decay phases (black) are shown in the second column. Temporal profiles 
of the plasma temperature (black) and the emission measure (gray) are shown in the third 
column. Evolution in the logT - log V EM plane (gray) with a fit to the flare decay slope 
((, black line) is shown in the fourth column. 



